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In this paper, we revisit one-channel time reversal (TR) experiments through multiple scattering 
media in the framework of the self-consistent multiple scattering theory. The hyper resolution and 
the self-averaging property are retrieved. The developed formalism leads to a deeper understanding 
of the role of the ladder and most-crossed diagrams in a TR experiment, and also establishes the 
link between TR and Coherent Backscattering (CBS). Especially, we show that when the initial 
source and the time reversal point are at the same location, the time-reversed amplitude is twice 
higher. Surprisingly, this enhancement is due to the ladder diagrams and not to the most-crossed 
ones, contrary to CBS. These theoretical predictions are confirmed by experimental results. The 
experiments are performed with ultrasonic waves propagating through a random collection of parallel 
steel rods. 

PACS numbers: 43.20-|-g,43.35.-|-d,73.20.Fz,42.25.Dd 



I. INTRODUCTION 

A few years ago, we have carried out the first experimental demonstration of the reversibility of acoustic waves 
propagating through a 2-dimensional random collection of scatterers, the ultrasonic equivalent of a Lorentz gas for 
particles T]. Typically, an ultrasonic source sends a broadband pulse into a medium consisting of thousands of parallel 
steel rods randomly distributed and immersed in water. Since the mean free path is much less than the sample 
thickness, strong multiple scattering occurs. The multiply scattered signals transmitted through the medium are 
recorded on a 128-transducer array, digitized and time-reversed. A part of the signal, the so-called "time reversal 
window" , is transmitted back into the medium. The time reversed wave is found to converge back to its source. 
Two aspects of this problem have been studied : the signal recreated at the source location (time compression) 
and the spatial focusing of the time-reversed wave around the source location. Dowling and Jacksonf^ have been 
pioneers in this field: three years before the first experimental observation, they predicted the peculiar property of 
super-resolution of time reversal in random media. But their approach was restricted to narrow-band signals and 
time-compression was thus not expected. Later we developed a phenomenological statistical model to describe both 
spatial and temporal focusing of a broadband pulse0,|3,0- Blomgren et al.0 applied the Green's function formalism 
to analyze the spatio-temporal focusing obtained by time reversal in a random medium with weak celerity fluctuations. 
However, their model was based on the parabolic (or paraxial) approximation(3| which assumes only weak scattering 
angles and thus ignores backscattering. In this latter paper, the reason why the focusing is observed for one realization 
of disorder is referred to as "pulse stabilization". It comes from the "self-averaging" behavior of the time reversal 
wave due to the broadband character of the initial pulse. 

In the present paper we revisit time reversal experiments in the framework of the multiple scattering theory. This 
theory describes statistical moments of the scattered field in transmission as well as in backscattering. It has been 
extensively applied in several fields of physics (optics, electronics, acoustics, ...). It is a general theoretical framework 
that successfully describes various multiple scattering phenomena: diffusive field fluctuations and their application 
to Diffusive Wave Spectroscopy (DWS)j3, 0, weak and strong localization ^3 short and long range intensity 
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FIG. 1: The two steps of a one-channel TR experiment through a muhiple scattering medium in a backscattering configuration, 
(a) forward propagation step, (b) backward propagation step. 



correlations (Ci, C2, C^)^^, etc... It was only very recently that arguments based on multiple scattering theory have 
been proposed in order to explain some results of TR in random mediafl^ . For simplicity, we focus in this paper 
on time reversal performed with one transducer. The average time-reversed field and its variance are worked out. 
Comparing the mean value of the TR field to its standard deviation leads us to understand the hyper-resolution as 
well as the self- averaging property of TR in random media. In addition to this, one of the most striking results of 
this paper is the following: the time reversal amplitude is predicted to be twice larger when the time reversal element 
and the initial source positions are the same. This effect is closely related to the well-known Coherent Backscattering 
Effect (CBS) ill m. Since the TR element is a receiver in the first step and a source in the second step of the TR 
process, the enhancement in TR focusing is due to the Ladder diagrams and not to the most-crossed ones, contrary 
to CBS! 

We first present some experimental results in order to introduce the questions that we address in this paper before 
presenting the theory in detail. The theoretical developments are based on the diagrammatic theory of multiple 
scattering. We apply it to determine the average amplitude and the variance of the TR field in both transmission 
and backscattering. In particular, we show that time reversal focusing depends on the so-called Vertex, which is the 
four-entry correlation function of the random wave-field. The study is limited to media whose Thouless number is 
much larger than one. In other words, the crossing paths inside the multiple scattering medium are neglected. The 
entire theory is explicitly developed in the time domain. We show that TR can be completely analyzed in terms 
of most-crossed and ladder diagrams. This general approach facilitates the study of two fundamental aspects of 
one-channel TR in multiple scattering media: hyper-resolution and self-averaging. We end the theoretical part by 
establishing a formal link between TR focusing and coherent backscattering. Finally, we show experimental results 
that confirm the theoretical predictions. 



II. A FIRST EXPERIMENT 



An acoustic wavefield can be time-reversed using a Time Reversal Mirror (TRM). A TRM is usually made of 
128 independent time reversal channels. Each channel consists of a piezoelectric transducer plugged to a digital 
emitting/recording electronics. In this article, we only deal with one-channel TR experiments performed through 
multiple scattering media. Indeed, the physics in such a configuration is already sufficiently new, rich and complex 
to devote the complete study to this topic. The schematic view of the one-channel TR setup is presented in Fig. Q 
The experiments have been carried out in a 2D-like multiple scattering medium made of thousands of steel rods. The 
sample is 35 mm thick and 300 mm large. There are 18.75 rods per square centimeter and the diameter of each rod is 
0.8 mm. The elastic mean free path is about 4mm [l^. The time reversal experiment is divided in two steps. First, 
a small emitter (here a piezoelectric transducer) located at point S generates a short ultrasonic pulse at a central 
frequency of 3.5 MHz (0.43mm wavelength) with a 100-percent bandwidth toward the multiple scattering medium. 
Another piezoelectric transducer acting as a small microphone records the time-dependence of the backscattered wave 
field, ^'(t), at point R. A typical backscattered signal is plotted in Fig. |21 Since the sample thickness is much larger 
than the mean free path, the huge time spreading (more than 100/is) of the initial pulse is due to multiple scattering. 
Secondly, the recorded signal is time-reversed and sent back into the medium by the transducer at point R now acting 
as a small loudspeaker. Finally a small transducer plugged to a waveform digitizer records the time-dependence of 
the back-propagated pressure field at point K. This transducer is fixed on a stepping motor in order to build the 
spatial map of the TR field around the initial source location. As expected, when K = S, a short pulse emerges 
that corresponds to the time-reversed initial pulse (cf. Fig. O and 0} This result is obtained for a single realization 
of disorder. In order to get an estimator for the average TR field, the experiment can be repeated for many other 
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FIG. 2: A typical time-evolution for the pressure field recorded at point R after the emission of a short l/xs pulse (at a central 
frequency of 3.5 MHz) at point S. 
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FIG. 3: The continuous line corresponds to the time compression recorded at the initial source position (K = S) for one 
realization of disorder. The dotted line is the average value over 100 realizations of disorder. The signals are normalized by 
the maximum of the average signal which occurs at t = 0. 




FIG. 4: Spatio-temporal focusing on a single realization of disorder (a) and averaged over 100 realizations (b). The vertical axis 
corresponds to the distance between the recording point (K) and the initial source position (S). The signals are normalized by 
the value of the averaged signal aX t = and K = S 
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FIG. 5: Average time compression on the initial source (K — S). Dotted line: the TR device location (R) is far from the 
initial source point (S), its amplitude is normalized to 1; continuous line: they are at same position (R — S). 
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FIG. 6: Full line: averaged time reversal focusing amplitude (at t = 0) and K = S versus distance between the initial source 
and the TR device. Dotted line: average backscattered intensity versus distance between the source and the TR device. 



configurations of the scatterer positions. From a practical point of view, averaging is achieved by translating the 
disordered medium. After averaging, the sidelobe level is decreased but the shape of the compressed pulse seems 
unchanged. This result illustrates the self-averaging property of broadband time reversal in random media. 

So far the source position and the time reversal channel were assumed to be far away from each other. Do the TR 
focusing characteristics change when points S and R coincide? An experimental answer comes from Fig. \S\ When 
S = K = R, the average TR field is almost twice as large compared to the set-up in which R and K = S are far 
away from each other. Moreover, if the average TR focusing amplitude is measured as a function of distance between 
points S = K and R, the resulting plot looks like the well-known Coherent Backscattering (CBS) peak (see Fig. (S)). 

The CBS is a well-known effect occurring in multiple scattering media, intimately related to the reciprocity property. 
When a source located at point S emits a short pulse towards a multiple scattering medium, the average backscattered 
intensity (i.e., J {^"ii (t)^^ dt where (•) represents the average over realizations of disorder) received at the source is 
twice higher than the one received far away. This is due to the constructive interference between a each path and its 
reciprocal counterpart, which can only occur at the source. This phenomenon has been observed in many different 
areas: in optics, in microwaves, in acoustics, etc... It was first predicted by Watson . d e Wolf Barabanenkov 
[l9|. Ten years later, the first experimental evidence of this phenomenon was reported 

Our preliminary results on time reversal focusing (Figs. [SlandEl) indicate that there is strong link between TR 
and CBS. We will later develop a full theoretical analysis to support this statement. Nevertheless, from now on, we 
propose an intuitive way to explain the link between CBS and TR based on simple physical arguments. Basically, 
when the source at S transmits a Dirac pulse, R records the Green's function of the medium between point S and R, 
G(S — > R;t). This field is time-reversed and G(S — > R, ~t) is sent back into the medium from point R to the initial 
source point, S. The back-propagated field on S is G(S R; —t) (g) G'(R S; t). The reciprocity symmetry of the 
propagation medium implies that the emission and receiving points can be exchanged in the Green's function, i.e., 
G(S R; t) = G(R S; t). Therefore the TR field is expressed as JG^{S^ R, r)dr at the focusing time {t = 0). 
This term is just the integrated intensity of the scattered field received at R for a source at S. Hence, on average, 
the dependence of the TR amplitude as a function of the TR device position is nothing else than the backscattered 
intensity pattern, i.e., the Coherent Backscattering enhancement. These preliminary and qualitative results will be 
supported rigorously in the rest of the paper. Theoretically we want to determine the average value and the variance 
of the time reversal field. The theoretical steps are the following: 
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FIG. 7: The two steps of a one-channel TR experiment in transmission through a multiple scattering medium, (a) forward 
propagation step, (b) backward propagation step. 




FIG. 8: The two steps of a one-channel TR experiment in a backscattering configuration through a multiple scattering medium, 
(a) forward propagation step, (b) backward propagation step. 



• Time Reversal in a complex medium on a single realization of disorder. 

— Some basics about Green's function formalism in disordered media are reminded 

— The TR field is deduced from the Green's function formalism. 

— Its square is calculated 

• Averaging over realizations of disorder: 

— General expression for the average TR amplitude in terms of the Vertex function. 

— General expression of the variance of the TR amplitude in terms of the Vertex function. 

• Approximation of the Vertex: introduction of the Ladder and Crossed diagrams contributions. 

— Analysis of TR in the transmission configuration (cf. Fig. [TJ. 

— Analysis of TR in the backscattering configuration (cf. Fig. (SI . 

• Discussion about the self-averaging property of TR in a multiple scattering medium. 

• The link between CBS and TR enhancement is established. 

The results will be confirmed by experiments presented in the last part of the paper. 

III. GENERAL THEORY OF TIME REVERSAL IN DISORDERED MEDIA 

A. Expression of the Green's function inside a multiple scattering medium 

In a medium filled with scatterers characterized by a scattering potential e, the Green's function G(ri„ — )■ Tout',^) 
between two points ri„ and rout at frequency uj is given bvpcl: 

G{rin ^ rouuuj) = Go(r.„, -> ro„t; cj) (1) 
+ J Goirin r;t)e(r)G(r ro„t;tj)rfr 
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Go is the Green's function without scatterers. If this expression is developed, it leads to the well-known Born 
expansion: 

G(ri„ rout;uj) = Go(rm -> Touu^j) + J Goinn -> r;<)e(r)Go(r ro„t;w)dr 

+ J Go(r„, ^ r; ^)e(r)Go(r ^ r'; ^)e(r')Go(r' ^ r^ut; '^)dvdv' + • • • (2) 

Each term of this development corresponds to an order of scattering: no scattering, single scattering, double scattering, 
... The aim of the multiple scattering theory is to find the average Green's function and its higher moments ((G), 
(GG*), etc..) from the previous expression. One can prefer to introduce an operator O defined such that|2lj|: 

G{Vin ^ Vout]^) — Ge(r,;„ ^ Tout; Cj) (3) 

+ y"Ge(r„ ri;tj)e(ri -> r2;w)Ge(r2 ro„f ; w)dr2dri 

where Ge is the average Green's function, i.e., (G). Ge is also called the Green's function of the effective medium. 
This definition implies that (8) = 0. Nevertheless the next moments of 0, (e.g. its autocorrelation) are different from 
zero. Experimentally in acoustics one has access to the time- dependence of the fields rather than their frequency 
dependence. So we prefer to work directly in the time domain. An inverse Fourier transform of Eq. yields an 
expression in the time domain: 

G(ri„ — > Tout; = Ge (fin ^ rout; t) (4) 

+ / Ge{Vin -> Yi]t)®Q{vi -> V2\t)®Ge{Y2 ^ Youut)dY2ri 



where ® is the convolution operator for time variables. 



B. Time reversed field 



In this part, we formalize the one-channel TR process in terms of the Green's functions. As seen in the preliminary 
experiment, the first step of a TR experiment begins by the emission of a short pulse, /(t), around time < = by the 
source localized at point S. The recorded field at point R is expressed as: 



*(t) ^G(S^R;t)(g)/(t) 

From Eq. jSJl, we have: 

*(i) = Ge(S^R;i)®/(t) 

+ j Ge(S -^Yi]t)® e(ri Y2]t) ® Ge(r2 R; t)dYidY2 



(5) 



(6) 



In the following, the first term on the right hand side of Eq. will be neglected. This is fully justified in a 
backscattering configuration: the effective Green's function has no specular reflection at the interface because the 
scatterers density is weak "2^ and the ambient fluid is the same in the scattering region and outside. In transmission, 
this term can also be neglected when the medium is thicker than several elastic mean free paths because the effective 
Green's function decays exponentially with depth in the multiple scattering medium. In other words, in thick multiple 
scattering media, almost all the incoming wave is scattered at least once while traveling across the medium. Without 
the first term on the right hand side in Eq. ©, the propagation from S to R can be seen as a three-step process: 
first the propagation in the effective homogeneous medium from S to ri , secondly the propagation inside the multiple 
scattering medium between ri and Y2 and finally the propagation from Y2 to R, along all possible paths from ri to 
Y2 within the scattering region. 

Next, we select a part of the scattered signal, referred to as the time- reversal window (TRW). This equals 



e{t) = A^[t)W 



t-tp 
At 



(7) 



W is the rectangle function [23] . A is an amplitude factor that takes into account a possible amplification, but can be 
set to unity without loss of generality. The TRW is centered at time to and has a duration AT. The selected signal 
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e{t) is time- reversed and transmitted back by point R through the medium and the field '^rt is recorded at point K 
(see Fig.|7|D and Eh): 



Combining Eq. lO, Eq. Q and Eq. (jHJl, we obtain 



(8) 



^'ijrW = J 7(1-1, r2,r3,r4;0'^/(-0 



W 



-t-tp 
AT 



e(ri ^ r2;-t) 



dridr2dr3dr4 



(9) 



X(i) 



where 



7(ri:i"2,r3,r4;0 = Ge(S -> ri; -t) Ge(R ^ rg; t) 

«) Ge(r2 ^ R; -t) (g) Ge(r4 ^ K; t) 



(10) 



To obtain Eq. ^ we have assumed that AT is larger than the duration of 7. 7 contains all the effective Green's 
functions and x takes into account propagation inside the disordered medium and the choice of the TRW. Expression 
(O for the time-reversed field in terms of Ge and the O matrix will serve us as a starting point to work out the average 
and the variance of the TR field. 



C. Average TR wavefield 

In Eq. lO, the randomness is contained only in x{t)- Hence, (x(0) is the key to (^_rt)- In Appendix IXI it is shown 
that: 



(X(^)) = 



y GT(ri,r2;r3,r4;T)M^ 



AT 



dr 



m 



(11) 



where 



GT(ri,r2;r3,r4;T) = (e(ri r2;T)e(r3 r4;T)) 



(12) 



To obtain Eq. (dTJ we assumed that in the frequency domain the correlation properties of the 0-operator only depend 
on the frequency difference. In the multiple scattering theory Ct is also referred to as the Intensity "Vertex" function 
that links the entry positions ri, ra to the output ones r2, r^. 

We conclude that the average TR field is intimately linked to the Intensity Vertex according to 



(*flT(0> = / 7(1-1, 1-2, rs, r4; i) ® /(-^) X 



J CT{ri,r2;r3,r4;T)W 



T-tp 



dr 



AT 
dridr2dr3dri 



Upon introducing the f function 



C.(Xi,X2,X3,X4;t) = yGe(Xi ^ri;-0®Ge(r2 ^X2;-i) 



Eq. lO reads: 



(g) Ge(X3 ^ ra; t) (g, Ge(r4 X4; t) (g, f{~t) 
X Ct (ri , r2 , ra , r4; T)dr 1 dr2dradr4 



(*fl.T(i)) = / er(S,R,R,K;i)W^(^^)dT 



(13) 



(14) 



(15) 



In the following, we will see that the ^ function plays a central role in the multiple scattering theory applied to TR. 
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From a practical point of view, we often deal with a short Time Reversal Window. When its duration AT is 
sufficiently small so that C-r (r i, r3, r4; r) is almost constant within this interval, the average TR field can be 
expressed as: 



AT6o(S,R,R,K,t) 



(16) 



AT^O 

The other interesting limit is when the Time Reversal Window contains the whole recorded signal ^'(t), in which case 



AT- 



(17) 



So far the average TR field has been worked out. To this end, we have assumed that the first term of the right hand 
side of Eq. ® is negligible. If we do not neglect this term, the term G'e(R ^ K; t) (g) Ge{S — > R; —t) (E) f{—t) must be 
added to Eq. Q as well as in the subsequent equations. To determine the variance of the TR field, we have to average 
the squared time-reversed field. The coherent terms are neglected under the same conditions as those explained in 
the previous section. Therefore the average squared field is: 



j [7(ri, r2, rg, r4; t - r) ® /(r - t)] W^.v'^, r^, r^; t-r')® /{t' - t)] 

(e(r3 - 
X (e(r3. 



(e(r3^r4;T)® [W^(-^ 
W{- 



r4' ; T' 



^^^)e(ri 

r'-tp ^ 



e(ri, 



-)]) 

-r2.;-r') 



X2{r,r') 

xdvidv2dv^dTidTi'dv2'dv^idvi>dTdT' 



(18) 



All effects induced by multiple scattering are governed by the function X2- Obviously the directly depends on 
the eight-point function in the frequency domain, (0(r3 r4;a;)9*(ri r2;cLi)9(r3' r4'; cj')9*(ri' T2'\ijj'))- In 
the following, we assume that the Thouless number of our medium is much larger than one so that all the random 
fields are jointly Gaussian and the cumulant theorem can be applied. Thus the eight-field correlation function can be 
expressed in terms of the four-field correlations. This approximation means that the probability for having crossing 
paths in the medium is weak[l^. The full calculation presented in Appendix IXI is tedious but without particular 
difficulty. Finally it comes: 



+ 6'(S,R,R,K,t)W^ 



e-t'(S,R,S,R,r)W^ 



2 / ^ *0 









t'=t T=2t 



AT 

e-t'(S,R,R,K,T)W^ 



6,(K,R,K,R, T)dT 



-t' - to 
AT 



(19) 



The subscripted (g operator notation is introduced in order to remove any possible ambiguity. f{t') ®t'=t 9{t') 
f^oo /(i')5(* ~ t')dt'. In the short and long TR window hmits respectively, Eq. ((T^ simplifies to: 



{^HTm - {^BT{t)}\^_^^AT J eto(S,R,S,R,r)6o+t(K,R,K,R,T)dT 



+Ar6o(S,R,R,K,T) <g> 6o(S,R,R,K,r) 

T=2t 



(20) 



i^RTit)') - {^RT{t))^) 



AT- 



C-t' (S, R, S, R, T)^t' (K, R, K, R, r)drdt' 
/"Ct'(S,R,R,K,T) «) £_^t'{S,R,R,K,T)dt' 

J T=2t 



(21) 



D. Approximation of the Vertex 

The Vertex can be split into 2 distinct contributions 1 1] : the irreducible one (U) and the reducible one (R). 

(e(ri ^ ra; t)e(ri, ^ ra- ; t)) = C/(ri, ra, ri- , rz- ; t) + i?(ri, ra, ri- , ra- ; t) (22) 
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FIG. 9: Diagrammatic representation of the expansion of the Vertex: (S) represents the single scattering contribution, (R) the 
"Ladder diagrams" and (C) the most-crossed diagrams that restore reciprocity 



The "Boltzmann approximation" consists in replacing (U) by the scattering by one scatterer (here noted S). This low 
density approximation turns (R) into an incoherent multiple scattering series called 'the Ladder diagrams" which obeys 
the radiative transfer equation. (R) itself does not obey reciprocity. Therefore the Boltzmann approximation fails in 
describing any reciprocity-dependent effect, such as Coherent Backscattering. Beyond the Boltzmann approximation, 
reciprocity is restored by adding "most-crossed diagrams" (C), which are irreducible, to the diagrammatic expansion 
(see Fig.EJ. Obviously, due to the straightforward link between the Vertex function and the £, function, (see Eq. 1^3), 
the same decomposition appears for ^: 

(Xi, X2, X3, X4, i) — ^^(Xi, X2, X3, X4, t) 

(Xi, X2, X3, X4, t) + (Xi, X2, X3, X4, t) (23) 

In the following, we shall neglect the single scattering contribution. The reciprocity symmetry means that all Green's 
functions are invariant upon permuting source and receiver positions : (G(ri ~^ r2',t) — G{t2 — > ri; t)). Consequently, 
in reciprocal media, a strong relation exists between the Ladder and the most-crossed diagrams; i?(ri, r2, r2' , ri' ; t) = 
C(ri, r2, Yy , Y2'\ t). The same relation holds for ^'^ and 

(Xi, X2, X3, X4, t) = (Xi, X2, X4, X3, t) (24) 

The solution for the total Vertex Ct (ri, r2; rs, r4; r) can be written F(ri, r2; t)5{Yi ~y^)5{y2 — r4). Under the diffusion 
approximation, F{Yi,Y2\t) obeys the diffusion eauation[2^. Using Eq. lfn|) gives 

e^(Xi,X2,X3,X4,<) = J Ge(Xl^Yi;-t)(g>Ge{Y2-^:^2;~t) 

(E) Ge(X3 ^Yi;t)(g) Ge{Y2 ^ X4; t) (g) f {^t)F {yi, Y2; to)dYidY2 (25) 

. This expression for ^ is valid in the diffusion approximation and takes into account the reciprocity. The next 
two sections of the article will be devoted to elucidating the respective role of the reducible and the most-crossed 
components of the vertex function in TR focusing. First we will consider the transmission configuration and next the 
backscattering configuration which is more complicated. 



IV. TRANSMISSION CONFIGURATION 



1. Average TR field 

In transmission, the points S and R are on opposite sides of the multiple scattering medium. The consequence 
is that ^^(S, R, R, K, t) is null. Indeed, the effective Green's functions are appreciable only on a skin layer whose 
thickness is one or a few mean free paths. If ri is in the skin layer at the sample input, then the effective Green 
function G'e(R Yi;t) is non zero but in this case Ge(S ri; —t) is negligible. Reciprocally, if Ge(R Yi;t) is non 
zero, Ge(S ri; —t) is zero. Therefore the average TR field (Eq. (|15() ') becomes: 

i^RTit)) = I (S, R, R, K; t)W{^)dr (26) 

That is to say that in transmission TR focusing properties come only from the existence of the most-crossed diagrams. 
A visual interpretation of the effect of averaging over realizations on TR propagation is proposed in Fig.^] The same 
kinds of conventions as used by Akkermans and Montambaux in ref j24l | are adopted here. For clarity only propagation 
inside the multiple scattering medium is represented. An acoustic path starts from ri and goes out at point Y2 with an 
associated. The wave is time-reversed. The new path starts from r3 and exits at Y4 with another phase. Time parity 
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FIG. 10: Paths that contribute to the TR mean field in transmission configuration. Dotted line indicates conjugated path. 




I II lU 

FIG. 11: Paths that contribute to the average squared TR field. 



transformation is formally equivalent to replace the forward path by its conjugate (see Fig. [Tn|l . The average process 
cancels out most of these paths. The only remaining paths are those that propagate along the same trajectory during 
the forward and the backward steps, (i.e. ri = r4 and T2 —v^, see Fig. I10|l . . Moreover the path and its conjugate 
have opposite propagation direction. It means that these paths will constructively interfere only if one of them can 
be replaced by its reciprocal. Thus the focusing only originates from the most-crossed diagrams that result from the 
reciprocity property. As the pairs of paths are linked by a convolution operator, they interfere constructively only at 
time t = Q. Usually in this kind of representation, a simple product operator replaces the convolution operator since 
one works in the frequency domain. Then constructive interferences occur at any time. 

Therefore in transmission and on average, the TR focusing only succeeds thanks to reciprocity. In other words, only 
the most crossed-diagrams contribute to the mean TR field. Usually the most-crossed diagrams are thought to be 
significant only in backscattering; here we have shown that they play a fundamental role in transmission as well when 
a time reversal experiment is performed. Now let us consider the statistical fluctuations around the mean value. 



2. Variance of the TR field 

Using the reciprocity principle (Eq. H24|l ) and the observation that in transmission ^^(S, R, R, K, vanishes, one 
infers that: 

(S,R,K,R,i) = (27) 

Thus the second term of the variance in Eq. Ijiyi) only originates from the Ladder contribution while the third term 
comes from the most-crossed diagrams. Again, a graphical interpretation is proposed in Fig. 1111 In the symbolic 
equation represented in Fig. there are now four paths in each term: two correspond (continuous line) to the 
forth and back propagation in the TR process and the other two symbolize the conjugate paths required to compute 
the square of the TR field. The three terms on the right hand side of the figure correspond to the 3 terms on the 
right hand side of Eq. l(TT))) . The first one denotes the square of the average value of the TR field. The second one 
corresponds to the contribution of the forward path (from ri to r2) with its conjugate together with the backward 
(from r3 to r4) path and its conjugate. The propagation directions of paths in each pair are the same. This means 
that the second contribution does not involve reciprocity of the medium. These two paths are linked together by a 
simple product operator, and is therefore a background contribution that slowly changes with time. The third term 
due to reciprocity only appears at the focusing (t « 0). 



V. BACKSCATTERING CONFIGURATION 



In the backscattering configuration, ^^(S, R, R, K, i) can be non zero when S and K are both close to R. In 
order to illustrate the theoretical results obtained in the backscattering configuration and to compare them later to 
experiment, we will consider an explicit solution for the ^ function. In the general case, the computation of ^ is very 
complex. However, it can be greatly simplified if we assume that all active elements (S, R and K) are in the far field 
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zone of the multiple scattering medium and that the initial pulse, f{t) has a bandwidth much smaller than its central 
frequency. Under these conditions, f is written as: 

e.(Xi,X2,X3,X4;t) = [F'(X4-Xi,X3-X2;r) 

+F'(X3-Xi,X4-X2;T)]/(-t) (28) 

where F' is the 2D-spatial Fourier transform of F' which is a function directly linked to F f Appendix Ib|) introduced in 
Eq. (|25l) as a solution of the diffusion equation with boundary conditions and calculated recently ,25] . The computation 
is quite long and complex and is out of the scope of this article. We have applied these results using parameters that 
correspond to our experimental ultrasonic experiment performed in a water tank. The central frequency of the initial 
pulse is 3.5 MHz, the sound speed in the homogeneous medium is c = 1.5 mm/ fis, the thickness of the multiple 
scattering medium is 35 mm. The diffusion coefficient is D ^ 3.2 mm^ / s, the elastic mean free path is ^ = 4mm and 
the transport mean free path is £* = 4.8 mmfT^. 



3. Average TR field 

From Eq. (|15|l and Eq. (|23|l show that the average time reversed field involves now both the ladder and the most- 
crossed diagram contributions. They are represented in terms of paths in Fig. 1121 The contribution / due to the 




/ // 

FIG. 12: Paths contributing to the average TR field in backscattering 



crossed diagrams corresponds to the one previously considered in transmission. However when the four points are 
at the same location, the second contribution (//) adds up to the previous one. It describes constructive interferences 
between a path and its conjugate that follows the same sequence of scattering events in the same order. As the 
direction of propagation is the same along both paths, this contribution would survive even if the reciprocity was 
broken. Both contributions show up only around time t = 0. The spatial focusing including these two terms is plotted 
in Fig. 1131 The ladder diagrams (contribution II in Fig. I12f) only contribute when S, R and K are sufficiently close 
from each other. The reciprocity principle links each ladder diagram to a crossed one. Therefore the amplitude of 
the TR peak is twice higher when the TR channel and the source coincide (see Fig. 114(1 . This behavior is very similar 
to the Coherent Backscattering Effect(CBS). More details are given in section IVIIII A surprising conclusion can be 
already drawn : when the reciprocity property of the propagation medium is broken, for instance by a strong flow. 




-a □ □ 



-5 □ -3 



FIG. 13: Numerical computation of the spatial dependence of the average TR wavefield at time t = G versus the distance 
between the source and the measurement points (vertical scale) and the distance between TR channel and the initial source 
(horizontal scale). The representations are in linear gray scale. The vertical curve is the plot of the field for R = S and the 
horizontal one for K = S. (a) and (b) correspond to TR focusing for time reversal windows respectively centered at 15 /is and 
55 jj,s. 
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FIG. 14: Evolutions of the average time-reversed peak amplitude recorded at the source (K = S) with respect to the center 
of the TR window (time to). Continuous line: TR channel is at the same position as the source (S = R = K). Dashed line: S 
and R are far away from each other. 



TR focusing still works as long as the initial source and the TR channel are at the same position. This is contrary to 
the intuitive idea that time reversal requires at least reciprocity of the medium. 



4- Variance of the TR field 

The expression of the variance of ^'_rt is quite complex. Indeed, the variance is composed of not less than 8 
terms. Nevertheless in the far field limit and for short TR windows, its expression is greatly simplified (Eq. H20|l and 
Eq. (EHIl): 

(*i?T W') - (*fl,T(i))' A^-^o^^'^^' °' + ^'(R - S, R - S; to)) 



x{F'{0,0;to+t)+F'{R~K,-R-K;to+t))AT J f{T)'^dT 



F'(K - S, 0; to) + - S, K - R; t^)] ' AT [f{-t) (g> f{-t)] (2 t) (29) 

The variance reduces to the sum of two terms. The first one represents a slowly varying contribution. If the source 
and the TR channel positions are identical (R = S) this term is twice higher than if they would have been far away 
from each other. The same conclusion applies when the TR channel and the recording point are identical (R = K). 
Thus when all three points coincide, we expect the variance to be enhanced by a factor four. 

The second term oscillates twice faster than the central frequency of the initial pulse. This term is significant only 
at the focus location (S) and at the focus time {t = 0). Moreover when the three points are all at the same location, 
this term is also enhanced by a factor four. Altogether the variance can be enhanced by a factor up to 8. In Fig. [T^ 
is plotted the spatial dependence of the variance at 3 different times for which the oscillating term is respectively 
minimum, zero and maximum. 

In conclusion, we have found three contributions to the variance of the TR field at backscattering. The first one is a 
spatially uniform background due to the ladder diagrams. In addition to this the most-crossed diagrams are responsible 
for a "diagonal enhancement" (R — K) and a "vertical enhancement" (R — S) (see Fig. I15|) . Now we have complete 
expressions and diagrammatic interpretations for the mean value and the variance of the TR field in transmission 
as well as in backscattering. In the next two sections, we apply these results to demonstrate "hyper-resolution' and 
"self-averaging" within the framework of the self-consistent diagrammatic approach. Finally we establish the link 
between TR focusing and Coherent Backscattering (CBS). 



VI. HYPER-FOCUSING 



One of the most striking effects in one-channel TR inside scattering or reverberant media is the hyper-focusing 
effect 0, 0, IS HE 113 ■ Resolutions of about one wavelength have been experimentally observed. In transmission, we 
can interpret this effect using Eqs. (|24|) . H25|l and H2t)|) . Indeed the average TR field can be seen as a converging field 
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FIG. 15: 3 Grey level representations of the variance. The vertical scale is the distance between the observer K and the 
source S. The horizontal scale is the distance between the TR channel R and the initial source S. The time reversal window 
is centered at to — 15 /xs. Figures (a), (b) and (c) represent the variance at 3 different times for which the oscillating term is 
respectively minimum, zero and maximum. 



emerging from an adaptive lens with an aperture function given by 0(ri) = /(/ F{^i,J^2',t)W[{t — to)/AT]dT) x 
Ge(R r2; t) (E) Ge{v2 — > R; —t) ® f{—t)dr2 where ri is the coordinate on the virtual lens. The maximum width of 
0(ri) is the transverse dimension of the sample. For an infinite slab, half-wavelength resolution is achieved with only 
a single channel TR! From a practical point of view, if the source is not really point-like, its angular spectrum sets 
the lower limit for the spatial resolution. In backscattering, the same formalism can be applied. 



VII. SELF-AVERAGING PROPERTY OF TR 



The time reversed signal at the focusing point (S) consists of a short pulse surrounded by side lobes that fluctuate 
from one realization of disorder to the other. We define the Signal-to-Noise ratio (SNR) at the focusing location 
(K=S) as the intensity of the peak (at t = 0) divided by the variance next to the peak. In other words, SNR = 
(^_RT(0))^/(5'flT(i)^)- Giving a general expression of the SNR is quite complex. However if the initial excitation 
function f(t) is narrow-band and frequency-dependent dissipation is negligible then Ge(S ~* ri; — (g) Ge(ri — > 
S; t) can be approximated by a Dirac S{t). In that case the SNR becomes for short and long time-reversal windows 
respectively: 

SNR = , (30) 

and 

<,NR = /^(O) {!F{r,,r2,T)drMrf 

AT^oo //2(t)dt /F(ri,r2,r)2dTdridr2 ^ ' 

Equation l|3U|) is equal to the ratio between the duration of the time reversal window and the one of the initial pulse 
f{t). As for Eq. H31|l . it represents the ratio of the time spreading of the initial pulse due to diffusion and the pulse 
duration. Hence for small TR windows the SNR ratio increases linearly with respect to the TR window length and 
reaches a plateau for large TR windows. This result, established here within multiple scattering theory, has been 
also found from phenomenological assumptions in several previous papers (e.g. Q). Expressions H3U|) and H31|) also 
show that the SNR increases as the pulse length decreases, i.e. when the bandwidth is enlarged. Broad band time 
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FIG. 16: Diagrammatic representation of the intensity (a) and the time-reversal amplitude at the focal spot (b). On the left 
part of the diagrams are represented the position of the 2 sources : twice S in (a), S and R in (b). On the right part of the 
diagrams are represented the positions of the 2 receivers : twice R in (a), R and S in (b). Only second order diagrams are 
drawn here. The left diagrams corresponds to Ladder contribution and the right ones to most-crossed contribution. 



reversal focusing in a complex medium tends to be a self-averaging process: the TR field obtained on one realization 
of disorder is a good estimator of its average field, provided that frequency-dependent dissipation is neglected. An 
experimental illustration has been provided in section ITU fFig. I^and Fig.^. 



VIII. TIME REVERSAL AT BACKSCATTERING VS. CBS 



To record the CBS in an acoustic experiment f |2f)ll28l |). a source located at S sends out a short pulse into a multiple 
scattering medium. The average reflected intensity is recorded at time to for each point R. 

/ = (G(S^R;to)') 

« 6o(S,R,S,R,t = 0) (32) 

This expression defines the spatial distribution of the baskscattcred intensity at a given time. Interestingly, the average 
TR field for S = K and at i = reads 



(*KT(t)> = AT6o(S,R,R,S,t = 0) (33) 

The expression of the average TR field is very close to the expression of the average backscattered intensity (Eq. I32|l . 
Yet, the exact equality is obtained only if the third and fourth entries are permuted. The physical signification of this 
permutation is fundamental. It implies that from TR to the CBS the role of the ladder and most-crossed diagrams 
are exchanged: in the CBS, the "background" and the "enhancement" intensity are respectively due to the ladder and 
the most-crossed diagrams, while in the average TR field, the background originates from the most-crossed diagrams 
and the enhancement from the ladders! 

This property appears more explicitly in the diagrammatic representations of the average intensity on the one hand 
and the average TR field amplitude on the other hand (Fig. I16|l . Let us recall the conventions of diagrammatic 
representations. First, the sources are on the left part of the diagrams while the receivers are on the right part. 
Roughly, a dot can be interpreted as a scattering point and two dots linked by a dashed line represent the same 
scatterer. A thick segment between two dots or between a dot and the one of the four outer positions symbolizes the 
effective Green's function; it is barred to indicate complex conjugation. The four outer positions are classified in two 
pairs : two positions form a pair if they are linked to the same scatterer by two Green's functions, one of them being 
conjugated. A rule exists to determine whether a diagram is significant or not. A diagram will be significant only if 
the two positions of each pair are sufficiently close and the contribution will reach its maximum when the two positions 
are identical. Hence, in the case of the intensity (Fig. 116b). the well-known result is retrieved : the ladders contribute 
whatever the distance between R and S whereas the most-crossed diagrams are only significant when R « S. In the 
case of TR, the role of diagrams are now exchanged (Fig. 116b): the crossed ones always contribute to the focusing 
amplitude whatever the distance between R and S. As to the ladders, they only arise when R w S! 

A direct consequence is that in a non-reciprocal medium the "enhancement" of the TR average field would remain, 
whereas the background would vanish. On the contrary, the enhancement of the backscattered intensity (i.e., CBS) 
would disappear, whereas the background would remain. 
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FIG. 17: Spatio-temporal representation of the TR field in gray scale. The vertical axis corresponds to the distance between 
the measurement point (K) and the initial source (S). The horizontal lower figure is the plot of the temporal focusing for 
K = S, and the vertical one is the spatial focusing at time t = 0. In both cases, the position of the Time Reversal Window is 
40/is and its duration is 4/i s. (a) is obtained when the TR channel (R) is far away from the initial source (S) and (b) when 
they are at the same position (R = S). The amplitudes are normalized such that the focusing amplitude in (a) equals 1. 




FIG. 18: Gray level representations of the average TR field at time t — versus the distance between the initial source and the 
measurement point (||S — K|| on the vertical scale) and the distance between the initial source and the TR channel (||S — R|| 
on the horizontal scale), (a) and (b) correspond to a 4 /xs long time reversal window centered at time to ~ 20 fis and to = 40^s 
respectively. 

IX. EXPERIMENTAL RESULTS 

The experimental setup was introduced already (see Fig. We now perform a "dynamic time reversal", i.e., we 
do not time reverse the whole signal but only a part of it defined by the Time Reversal Window W {{t — to) / AT) as in 
Eq. Q. The experimental results we present here have been obtained in backscattering, which is the most interesting 
configuration. First we focus on the average time reversed field, then on its variance. 

A. Average 

As in the preliminary part, the source (S) sends a 1-fis pulse at Z.bMHz central frequency. The backscattered 
ultrasonic field is recorded at point R. A short time reversal window centered at t^ = 40/Lts and AT = 4/Lts long 
is selected, time-reversed and transmitted back towards the medium. This experiment is repeated for one hundred 
realizations of disorder. We have recorded the average TR field for two extreme situations: either the one-channel 
TR device and the initial source are far away from each other fFig.llTb). or they are at the same location fFig. 117b). 
First we observe that the hyper-focusing property of TR through complex media is achieved: the focal spot is very 
thin [2mm i.e., 4 wavelengths) whereas no spatial focusing would occur in free space with one single TR channel. 

Furthermore we notice that the amplitude of the pulse is twice as large as when the TR channel is at the same 
position as the initial source (R = S). This is a manifestation of the coherent backscattering effect. As predicted by 
our theoretical study (see Fig. ll3|) . we observe in Fig. ^| that the TR focusing amplitude versus the distance between 
S and R is peaked like the CBS. The width of the peak narrows as to increases. In Fig.^|is plotted the dependence 
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FIG. 19: Amplitude of the average TR field (recorded at time t=0 and at the focal position, i.e., S — K) versus the position 
of the TR window (to). Continuous line (respect, dotted line) corresponds to the case where the TR channel is far from the 
source whereas the dotted line corresponds to the situation where the initial source point and the TR channel are identical. 
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FIG. 20: Average TR focusing peak amplitude at the source (S — K) versus the distance between the source and the TR 
channel (fr — rs) versus the position to of the TR window. At each time, the field is normalized such that the field equals 
1 when S and K are far apart. The top plot shows two snapshots of the bottom map at times to = 17^s (dashed line) and 
to = 57/is (dashed dotted line) 

of the peak amplitude with respect to the TR window position when S = R and when S 7^ R (Fig. I19|l . We observe 
the behavior predicted by the theory (see Fig. I14|l . Especially we observe that the factor 2 is measured whatever the 
Time Reversal Window as soon as there are at least two scattering events. Finally Fig. [201 presents a map of the peak 
amplitude with respect to the postion of the TR window (to) and the distance between the source (S) and the TR 
channel (R). Its width decreases as 1/y/Dto as to increases, like the dynamic CBSp^. 

B. Variance 

Figure El shows a map of the average TR field when the source point S and the single TR channel R are far 
away from each other. The corresponding map of the variance is shown in Fig. 1211 As predicted by the theory, two 
contributions add up to a flat variance that was normalized here to 1. Upon varying the position of the receiver 
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FIG. 21: Gray level representation of the variance of the TR field versus time t (horizontal scale) and versus the distance 
between the measurement point K and the initial source S (vertical scale). The TR channel is located at —4.6 mm from 
the source. The horizontal plot is the time evolution of the variance at the focus (S — K). The vertical plot is the spatial 
dependence of the field at t = 0. The TR window is AT = 2//s long and centered at time to = 40/xs. 
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FIG. 22: Variance of the TR field outside the focusing time. The vertical axis corresponds to the distance between the TR 
channel (R) and the source (S). The horizontal axis is the distance between the measurement point (K) and the source (S). 
The 2fis long TR windows are centered around time 20/is (a) and 40/is (b) respectively. 

K, the first contribution appears around the position of the TR channel R. This contribution is seen to be roughly 
constant in time. The second contribution only occurs at the source position S and for times around the focusing 
time t = 0. As was expected, this contribution oscillates twice as fast as the central frequency of the initial pulse. In 
both case, the maximum enhancement factor is 2. Finally the experimental results shown in Figs.|23and 1231 confirm 
the theoretical result that was plotted earlier in Fig. 

X. CONCLUSIONS 

One-channel time-reversal experiments have been revisited within the framework of the multiple scattering theory. 
On the one hand, previously established results such as hyper-resolution and self-averaging have been rigorously 
demonstrated within this theoretical framework. On the other hand, the link between one-channel TR and CBS has 
been established. Especially, we have shown that when the initial source and the time reversal point coincide, the 
time-reversed amplitude is twice as large. Surprisingly, this enhancement is due to the ladder diagrams and not to 
the most-crossed ones, contrary to CBS. These theoretical predictions have been confirmed by experimental results 
which have been obtained with ultrasonic waves propagating through a random 2D collection of parallel steel rods. 
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FIG. 23: Variance of the TR field at the focusing time. The vertical axis corresponds to the distance between the TR channel 
R and the source S. The horizontal axis is the distance between the measurement point K and the source S. The TR windows 
are centered around time 20yU,s (a) and 40/Lts (b) respectively. 

The generalization to multiple-channel TR and to random media subject to spatial long-range correlations will be the 
object of further studies. 
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APPENDIX A: MEAN VALUE OF x AND X2 



Let A and B be two time-dependent, real-valued, random fmictions that in the text correspond to 0(ri — )■ and 
Q{r^ r^) respectively. We want to find the expression for the correlation between these two functions from their 
frequency correlations. From the definition of the Fourier transform, the time correlation is 



{A{t)B{t')) = J {A{ij)B*{u;'))e'^-'^'+'^'''^du;du;' 



(Al) 



Here and later, the integrations are implicitly performed from —oo to oo. We have assumed that the correlation 
function in the frequency domain Ca,b depends only on the difference 5lo — lu — lo' . CA,Bit) is the inverse Fourier 
transform of Ca,b{S'-^)- To calculate the mean TR field, we need to compute 



{B{r)A{T-t))w C L ^° )dr 



AT 



ixit)) - 

Using Eq. (|Aip . the last expression reads 

ixit)) = 

Evaluation of the variance of the TR process implies to work out: 

(X2(r,r'))= / (BirM{ri-r)B'{r[)A'ir[~r')) 



CAM{r)W{^)dT 



m 



(A2) 



(A3) 



(A4) 



^ ^ AT AT ^ ^1 

Assuming A,B,A',B' are jointly Gaussian variables, this becomes 

{Bin)Ain-r)B'ir[)A'ir[~r'))^ 

{B{r,)A{ri-r)) {B' [r'M' {t[ - r' )) + 
{B{r,)B'{T[)){A{T,-r)A'{T[-r')) + 

CB,B'{-^l)CA,A'iri-r)5{Ti-T[)5{Ti-T[-T + T') 

{B{TM\r[ - r')> {B'ir'Miri - r)) (A5) 
^ V ' 

CB,A'{'^l)CB',A(-ri-'<-)S{Tl-T[+T')5{T[-Tl+T) 

Hence, (X2) is the sum of tree terms: 

(X2(t,t')) = Xi(r,r')+X2(T,T')+X3(r,r') 



where 



Xi(r,r') ^ (x(r)> (x(r')) 
X2(t,t') = (cbm'{t)® 



CA,A'i-^W{ 



AT 



yi-T + r') 



X3(r,T') = 



Cb'A~tW{ 



-T - to , 

AT ' 



S{t' 



(A6) 

(A7) 
(A8) 

(A9) 



These expressions can be greatly simplified if we assume that the correlation function is constant inside the duration 
of the Time Reversal Window AT 



X2(T,r') « ATCA,A'{t„)CB,B'{r + to)5i-T + t') 
X3(r,r') « A{^)CB,A'{to)CB'Ato)6{T' + t) 

where A is the triangle function, i.e., A(a;) = 1 — |x| if \x\ < 1 and if |a;| > 1. 



(AlO) 
(All) 
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APPENDIX B: FAR FIELD APPROXIMATION 



The far field approximation leads to write the convolution of two effective Green's functions as 

Ge(Xi^ri;-i)®Ge(X2^ri;t) = p(5(t-^i^^^^-^^)x I e ^when2>0 (bi) 

ca when z < 

where a is the distance between the plane that contains Xi, zi and xi which are respectively the normal and parallel 
components of the position vector with respect to the plane interface of the multiple scattering medium, p is a constant 
equal to 1/a* in 3D. £ is the elastic mean free path which is assumed to be independent of frequency (resonances are 
ignored), /ii — cos{9i) where tan(0i) = Xi/a. Then Eq. I|14() can be rewritten as 

6o(^l'^2,X3,X4;t) = 

S{t + xi— + X2— ^)F'(xi, X2; to)dxidx2 <^ f{-t) 

ca ca 
. ' 

c 

6{t + xi— + X2 -)F' {xi,X2;to)dxidx2 <^ f{~t) 



where F' is given by: 

F'{xi,X2;to) = p J e-''^^/^'-'e-'''^l'''F{xi,X2,zx,Z2-M)dzxdz2 (B2) 

The integrations are always implicitly performed from — cx3 to cxd. A narrow band pulse can be treated as almost 
monochromatic. Then ^ can simply be expressed in terms of the spatial Fourier transform of F' : 

Cto(^i'^2,X3,X4;i) = [F'(X4 — Xi.Xa — X2;to) 

c 

+ i?'(X3-Xi,X4-X2;<o)]/(-i) (B3) 

R 

where the 2-D Fourier transform is: 

[ e-'^^F'{xi,x2;to)e-'^^dxidx2 (B4) 
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